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Abstract 

The decay K~ — > TT°fi~v has been studied using in-flight decays detected with the 
"ISTRA+" spectrometer. About 540K events were collected for the analysis. The A + 
and Ao slope parameters of the decay form-factors f+{t), fo(t) have been measured : 
A + = 0.0277 ± 0.0013(stat) ±0.0009(syst), A = 0.0183 ± O.OOll(stat) ±0.0006(syst), 
and d\o/d\+ = —0.348. The limits on the possible tensor and scalar couplings have 
been derived: /t//+(0) = -0.0007 ± 0.0071, / s //+(0) = 0.0017 ± 0.0014. No visible 
non-linearity in the form-factors have been observed. 



1 Introduction 

The decay K — > / uz/7t (K m3 ) provides unique information about the dynamics of the strong 
interactions. It has been a testing ground for such theories as current algebra, PCAC, Chiral 
Perturbation Theory (ChPT). In this paper we present a high-statistics measurement (~ 537K 
events) of the Dalitz plot density in this decay. This study has a particular interest in view of 
new two-loop order (p e ) calculations for K; 3 in ChPT [1]. 

The K M 3 decay is also known to be a key one in hunting for phenomena beyond the Standard 
Model (SM). In particular, significant efforts have been invested into T- violation searches, by 
the measurements of the muon transverse polarization ot [2], as well as into searches for the 
non-SM contributions into the decay amplitude [3]. 

In our analysis we present new search for scalar (S) and tensor (T) interactions by fitting 
the K^ 3 Dalitz plot distribution, similar to the procedure used in the K e3 decay studies[4]. 

2 Experimental setup 

The experiment has been performed at the IHEP 70 GeV proton synchrotron U-70. The 
experimental setup "ISTRA+" (Fig.l) was described in some details in our paper [5]. 




Figure 1: Elevation view of the "ISTRA+" detector. 



The setup is located in a negative unseparated secondary beam. The beam momentum is 
~ 25 GeV with Ap/p ~ 1.5%. The admixture of K~ in the beam is ~ 3%. The beam intensity 
is ~ 3 • 10 6 per 1.9 sec. of the U-70 spill. The beam particles are deflected by the beam magnet 
Mi and are measured by BPC\ -j- BPC^ proportional chambers with 1 mm wire spacing. The 
kaon identification is performed by Co -j- C2 threshold (7-counters. 

The 9 meter long vacuumed decay volume is surrounded by 8 lead-glass rings LG\ 4- LGg 
which are used as the veto system for low energy photons. The photons radiated at large angles 
are detected by the lead-glass calorimeter SP 2 . 

The decay products are deflected by the spectrometer magnet M2 with a field integral of 
1 Tm. The track measurement is performed by 2-mm-step proportional chambers (PC1S-PC3), 
1-cm-cell drift chambers (DC\ -7- DC3), and by 2-cm-diameter drift tubes (DT\ 4- DT4). Wide 
aperture threshold Cerenkov counters (C3 and C4) are filled with helium and are not used in 
these measurements. 
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The photons are measured by the lead-glass calorimeter SPi which consists of 576 counters. 
The counter transverse size is 5.2 x 5.2 cm and the length is about 15 X . 

The scintillator-iron sampling hadron calorimeter HC is subdivided into 7 longitudinal sec- 
tions 7x7 cells each. The 11x11 cell scintillating hodoscope is used for the improvement of the 
time resolution of the tracking system. MuH is a 7x7 cell scintillating muon hodoscope. 

The trigger is provided by Si S 5 scintillation counters, C 4- C 2 Cerenkov counters, and 
the analog sum of amplitudes from last dinodes of the SP 1 : 

T — S± • S2 ■ S3 • S4 ■ Co • C\ ■ C2 • S$ ■ S(jS'Pi), 

where S4 is the scintillator counter with a hole to suppress the beam halo, S5 is the counter 
located downstream the setup at the beam focus. This part of the trigger is intended to identify 
beam kaons and to kill undecayed particles. It is designed on purpose, in a very simple way, to 
avoid any bias. E(S'Pi) requires that the analog sum of amplitudes from the SPi be larger than 
~700 MeV - the MIP signal. The last requirement serves to suppress the dominating K — > \iv 
decay. A part of events (10%) which do not satisfy the Y,{SP\) requirement is also recorded to 
provide the information for muon identification studies. 

3 Events selection 

During run in Winter 2001, 332M events were logged on tapes. This statistics is complemented 
by about 130M MC events generated with Geant3 [6] Monte Carlo program. The MC genera- 
tion includes a realistic description of the setup with decay volume entrance windows, tracking 
chambers windows, chambers gas mixtures, sense wires and cathode structures, Cerenkov coun- 
ters mirrors and gas, the shower generation in EM calorimeters, etc. 

The data processing starts with the beam particle reconstruction in BPC\ -4 BPC4. Then 
secondary tracks are looked for in the decay tracking system and events with one good negatively 
charged track are selected. The decay vertex is reconstructed by means of the unconstrained 
vertex fit of the beam and decay tracks. 

A clustering procedure is used to find showers in the SPi calorimeter, and the two- 
dimensional pattern of the shower is fitted with the MC-generated patterns to reconstruct 
its energy and position. 

The muon identification is done using the information from electromagnetic and hadronic 
calorimeters. First of all, the energy deposition in the SP\ associated with the track (counted 
in the 3x3 matrix around the track extrapolation to the SPi) is required to be less than 500 
MeV. This cut is intended to suppress the electron tracks. The sum of ADC counts from the 
HC counters associated with remaining tracks is demanded to be less than 200 (see Figure 2). 
And, finally, the ratio of the associated ADC signals in the last three layers of HC to the total 
associated ADC sum to be greater than 0.05 is required (Figure 3 shows this value for the 
tracks which pass the first two selection criteria). 
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Figure 2: The ADC sum in HC for the Figure 3: The ratio of the track-associated 
track-associated cells. ADC signals in the last three layers of HC 

to the total associated signal. 

The figures 2 and 3 are obtained with the clean 7r~ and \i~ samples selected from the real 
data. The pion data sample is composed from selected K~ — > 7t~tt° decays, and the muon one 
from the K~ — > i^v^ decays. 

The efficiency of the muon identification and the probability of the 7r — > fi assignment were 
found to be 88% and 0.03 respectively. 

The events with one charged track identified as muon and two additional showers in the 
SPi are selected for further processing. 

The selected events are required to pass 2C K — > fiu-K fit, with a probability of the fit 
P^ > 0.005. The angle between 7r° and yT in the kaon rest frame after 2C fit was found to 
be a good variable for the further background suppression (see Figure 4). The background 
from the surviving K~ — ► 7r~7r° events is concentrated at cos 9 ~ —1, and the selected cut 
cos^t,-^ > —0.95 removes practically all the background. The missing energy E v = Ek—E^—E^o 
after the angular cut is shown in Figure 5. The signal Monte-Carlo events for Figures 4 and 
5 are weighted with the matrix element where we use A + = 0.0286 (fixed from our K e % 
measurements [4]) and A = 0.017 (from the ChPT 0{p A ) calculations [7]). 

We estimate the surviving background contribution to be around 0.3%. 
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Figure 4: The cosine of the tt - p angle in Figure 5: The E v compared with MC. 
the kaon rest frame after 2C fit. The points 
with errors are data and the solid histogram 
is MC. The shaded area shows the 
background contribution. 



4 Analysis 

After the procedure described in the previous section, 537K events are selected in the real data. 
The distribution of the events over the Dalitz plot is shown in Figure 6. 

The most general Lorentz- invariant form of the matrix element for the K~ — ► Eutt decay 
is [8]: 

M= ^^n(p,)(l+7 5 )[2m K / s -[(P^ + P.) Q / + + (P^-P w ) a /„] 7 a +^^a^P^] W ( W ) (1) 

It consists of scalar, vector, and tensor terms. The f± form-factors are the functions of t = 
(P K — P T ) 2 . In the Standard Model (SM), the W-boson exchange leads to the pure vector term. 
The scalar and/or tensor terms which are "induced" by EW radiative corrections are negligibly 
small, i.e nonzero scalar or tensor form-factors would indicate the physics beyond the SM. 

The term in the vector part, proportional to /_, is reduced (using the Dirac equation) to 
the scalar form-factor. In the same way, the tensor term is reduced to a mixture of the scalar 
and vector form-factors. The redefined vector (V) and scalar (S) terms, and the corresponding 
Dalitz plot density in the kaon rest frame (p(E n ,Ei)) are [9]: 

p(E n ,Ei) ~ A ■ \V\ 2 + B ■ Re(V*S) + C ■ \S\ 2 (2) 
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V = f + + (mi/m K )f 7 



nij 2E, EA 



S = f s + (m l /2m K )U+ i + ^ r -_L--Jl f T 

A = m K (2E l E u -m K AE^)-mf(E v --AE 7r ) 
B = mim K (2E I/ - AE W ); E u = m K - Ei - E w 

max T7i . Tpmax _ m K ~ m l + m 7i 



C = m 2 K AE n , AE W = E™ ax - E n ; E w 



2rriK 



Following [7] the scalar form-factor fo is introduced: 

/o(t) = /+(*)+ 2 t 2 /-(*), 
and we assume, at most, the quadratic dependence of /+, f on t: 

/+(*) = / + (0) (l + X + t/ml + A^/m*) , /o(*) = /+(0) (l + U/ml + \' Q t 2 / 
Finally, one gets from Eq. (3): 

/- = / + (o)^s^ • f Ac - a + + -^(a; - a;)) 

ml \ mi ^ J 
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Figure 6: Dalitz plot for the selected K — > [lvtx events after the 2C fit. 



The procedure of the extraction of the form-factor parameters starts with the subdivision 
of the Dalitz plot region y = 0.425 0.955; z = 0.545 1.025 into 40 x 40 bins. 

The signal MC was generated with the constant matrix element and we have to calculate 
the amplitude-induced weights during the fit procedure. One can observe that the Dalitz-plot 
density function p(y, z) of (2) can be presented in the factorisable form, i.e 

p(y, z )= F a(K^'+^oA' JsjT)-K a (y,z), (6) 

a=l,18 

where F a are simple bilinear functions of the form-factor parameters and K Q (y,z) are the 
kinematic functions which are calculated from the MC-truth information. For each a, the sums 
of K a (y,z) over events are accumulated in the Dalitz plot bins (i,j) to which the MC events 
fall after the reconstruction. Finally, every bin in the Dalitz plot gets 18 weights W a (i,j) and 
the density function r(i,j) which enters into the fitting procedure is constructed: 

r(i,j)= ]T F a (X + ,X' + ,X ,X' J s J T ).W a (i,j) (7) 

a=l,18 

This method allows one to avoid the systematic errors due to the "migration" of the events 
over the Dalitz plot due to the finite experimental resolution and automatically takes into 
account the efficiency of the reconstruction and selection procedures. 

To take into account the finite number of MC events in the particular bin and strong 
variation of the real data events over the Dalitz plot, we minimize a —C function defined as 
[10]: 

-C = 2 n 3 ln 

3 

where the sum runs over all populated bins, and rij, rj and rrij are the number of data events, 
expected events and generated Monte Carlo events respectively. For large rrij Eq. (8) reduces 
to the more familiar expression 

j 

The minimization is performed by means of the "MINUIT" program [11]. The errors are 
calculated by "MINOS" procedure of "MINUIT" at the level AC = 1, corresponding to 68% 
coverage probability for 1 parameter. 

5 Results 

A fit of the fT M 3 data with f$ = fr = X' + = X' = gives the following result for A + and Ao: 
A + = 0.0277 ± 0.0013; A = 0.0183 ± 0.001. The A + - A correlation parameter is found to be 
dX /dX + = -0.348. The total number of bins is 1054 and xV ndf = i- 008 - The quality of the fit 
is illustrated in figures 7 and 8 where the projected variables y = 2E fl /mK and z = 2E 7T o/rriK 
are presented. 
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Figure 7: Y distribution. 
The points with errors are the real data 
and the shaded area - signal MC. 
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Figure 8: Z distribution. 
The points with errors are the real data 
and the shaded area - signal MC. 



The value A+ = 0.0277±0.0013 is in a good agreement with that extracted from the analysis 
of our K e3 data [4]: A+ = 0.0286 ± 0.00054(stat) ± O.OOOQ(syst) (the statistical error ±0.0008 
presented in [4] was obtained with AC = 2.3, corresponding to 68% coverage probability for 
2-parameter fit), i.e, our data do not contradict \i — e universality. 

In addition to the fits described above, Table 1 represents the fits with possible nonlinear 
terms in / + and fo (Eq. 4) as well as the fits with tensor and scalar contributions (Eq. 1). 

Every row of the Table 1 represents a particular fit where the parameters shown without 
errors are fixed. The second row shows a fit where the nonlinearity is allowed in f+{t). One can 
observe A + — A' + correlation that results in the significant A + errors enhancement and visible 
shift of A + and Ao parameters. The fitted value of A+ is compatible with zero, while we can 
not exclude some nonlinearity. The third row represents a fit with the value of A+ parameter 
extracted from the analysis of the data on pion scalar form-factors [1]: X' + = 3.2-m 4 = 0.001063. 
In a similar way A parameter is strongly correlated with Ao and is compatible with zero (row 
4). 

We do not see any tensor contribution in our data (row 5). The last row of the Table 1 
represents a search for the scalar contribution. As one can see from the Eq. (2), the fs term is 

2 2 

100% anti-correlated with V-A contribution (m^/2m K )f_, where /_ = / + (0)(A — A + ) my 
i.e an independent estimate of this term is necessary. A possible way consists in fixing Ao at 
the value calculated in the 0(p 4 ) ChPT: A$ h = 0.017 ± 0.004 [7]. The error ±0.004 in the 
theoretical prediction induces an additional error of ±0.0053 in f$/ f+(0)- 
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0.0277 ±0.0013 
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Table 1. The fits. 

Different sources of systematics are investigated. We allow variations of the muon selection 
cuts, angular cut and 2C-fit probability cut. The Dalitz plot binning, signal and background 
MC variations are also applied. 

The resulting systematic uncertainties are as follows: 

• AA+ = 0.0009 and AA = 0.0006; 

• A/ T //+(0) = 0.002 and A/ s //+(0) = 0.0009 

6 Summary and conclusions 

The K~ 3 decay has been studied using in-flight decays of 25 GeV K~ detected by the "ISTRA±" 
magnetic spectrometer. 

The A + parameter of the vector form-factor is measured to be: 

A+ = 0.0277 ± 0.0013 (stat) ± 0.0009 (syst) 

The Ao parameter of the scalar form-factor is defined: 

A = 0.0183 ± 0.0011 (stat) ± 0.0006 (syst) 
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The comparison of the A + parameter with that obtained from our K e3 data shows e — jj, 
universality. 

It is, at present, the best measurement of these parameters. It is in a reasonable agreement 
with 0(p 4 ) ChPT prediction as well as with recent Ao measurements from the T(K^)/T(K e3 ) 
ratio [12]. 

Possible quadratic contributions in the vector and scalar form-factors are compatible with 
zero, further studies are necessary to perform a detailed comparison of our data with 0(p 6 ) 
ChPT calculations [1]. 

The limits on possible tensor and scalar couplings are derived from the combined fit: 

/ T //+(0) = -0.0007 ± 0.0071 (stat) ± 0.002 (syst); 
f s /f+(0) = 0.0017 ± 0.0014 (stat) ± 0.0009 (syst) ± 0.0053 (theor) 

The work is supported by the RFBR grant N03-02- 16330. 
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